Viable Palatini-f(R) cosmologies with generalized dark matter 
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We study the formation of large-scale structure in universes dominated by dark matter and driven 
to accelerated expansion by f(R) gravity in the Palatini formalism. If the dark matter is cold, 
practically all of these models are ruled out because they fail to reproduce the observed matter 
power spectrum. We point out that if the assumption that dark matter is perfect and pressureless 
at all scales is relaxed, nontrivial alternatives to a cosmological constant become viable within this 
class of modified gravity models. 



I. INTRODUCTION 

The f(R) theories of gravity have been intensively ex- 
plored in recent years as possible alternatives to dark 
energy [lj. This class of extensions of general relativity 
is defined via a seemingly simple generalization of the 
Einstein-Hilbert action by allowing nonlinear interactions 
of the Ricci scalar as 



16nG 



-gf (g^R^(r)) + s n 



(1) 



where S m is the matter action. When T is taken to 
be the metric-compatible (Christoffel) connection, the 
action (T) represents the so called metric f(R) models 
@, H, UH, Hi- Though a realistic universe expansion 
history may be reconstructed from these models 0, @|, 
they seem to be ruled out as alternatives to dark en- 
ergy because of their consequences to the Solar system 
physics 

see however [l2| • If the connection T 
is promoted to an independent variable, the action (fTJ) 
represents the f (R) models in the so called Palatini for- 
malism [l3|, [Til, EE Gil E3, fl8|, w hich instead appear to 
pass the Solar system tests 19|. The viability of the 
small-scale limit of these Palatini-/ (R) models has been 
an issue of debate [20l l2lj ]. A recent claim of violations 
of the equivalence principle [22| . though partly erronous 
[23| . has resurfaced doubts about the physical tenability 
of the Palatini-/ (R) gravity [H, H3|- These have been al- 
ready considerably clarified by Kainulainen et al [Tl| , and 
the problems which perhaps remain at tiny scales could 
also dissipate if one entertains these models as a macro- 
scopic limit of a possibly more fundamental description 
of spacetime with the aid of truly metric-affine degrees 
of freedom 



25] 



This suggests it worthwhile to seriously reconsider the 
cosmology of the Palatini-/ (R) models. Expansionwise, 
they can generate a viable sequence of radiation dom- 
inated, matter dominated and accelerating era match- 
ing with the constraints, as has been shown for var- 
ious parameterizations of the function f(R), most of- 
ten with some power-law forms (with one, two or three 
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powers of R), but also with square-root, logarithmic 
and exp onential forms for the curvature correction terms 
[H, H3 IHj [Hj H3| ■ However, the inhomogeneous evolu- 
tion present in any realistic universe has not been suc- 
ccsfully reconciled with observations in these models. At 
the background level the additional derivative terms in 
the field equations due to nonlinearity in f(R) can play 
the role of an effective smooth dark energy. Meanwhile 
the spatial gradients of these extra terms cannot be ne- 
glected for cosmological perturbations [3l| which then 
assume unusual behaviour that is at odds with the ob- 
served distribution of galaxies [12] and with the cosmic 
microwave background spectrum (33j . even if the non- 
linear part of the action is exponentially suppressed at 
late times [34j |. In quantitative terms, parameterizing 
the f(R) — R ~ Rr, data analysis constrains \/3\ < 10~ 5 
[3lL [33] ■ Thus, while these models are allowed by local 
experiments, they are extremely tightly constrained by 
cosmological data. 

In this paper our purpose is to investigate under which 
conditions these cosmological constraints on the Palatini- 
f(R) could be loosened. Our approach is to allow gen- 
eralized dark matter (GDM) with possible (isotropic or 
anisotropic) pressures in the matter sector. The grav- 
ity sector we assume to be given by a general nonlinear 
function f(R) appearing in Eq.|T|), when understood in 
the Palatini formalism which yields second order field 
equations (in contrast to the "fourth order" metric for- 
malism). Note that a linear Lagrangian f(R) ~ R + 2 A 
corresponds to general relativity with a cosmological con- 
stant A. Note also that here we assume the nonlinear part 
of the function f(R) to drive the cosmic acceleration so 
that there is no need for dark energy. One could attempt 
to eliminate also the need for dark matter with nonlin- 
ear gravity 35]. That is not our approach here, but our 
results could be relevant for such pursuits too. 

The paper is organized as follows. We generalize the 
cosmological equations in section [TT] and apply them in 
section IIIII to construct an /(i?)GDM scenario, where 
anisotropically stressed dark matter has (up to a nor- 
malization) the ACDM (cold dark matter model with a 
cosmological constant) matter power spectrum, in a uni- 
verse with its the background expansion identical to the 
corresponding /(i?)CDM model. Therefore we can con- 
clude in section IIVI that with specific assumptions about 
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the properties of dark matter, the large-scale structure 
in these models can be consistent with observations. Ap- 
pendix |A] contains the full linearized evolution equation. 

II. COSMOLOGY WITH GENERALIZED DARK 
MATTER 

In the ACDM model the CDM component consists 
possibly of weakly interacting massive particles or some- 
thing else which should be at cosmological scales very 
accurately approximated as an exactly pressureless and 
perfect fluid. These properties are deduced from the ob- 
servations of large-scale clustering properties of matter, 
assuming the gravitational dynamics to be governed by 
GR. The so called hot dark matter scenario, as an ex- 
ample, is excluded because of the finite pressure in a hot 
component inhibits matter from clustering at subhorizon 
scales efficiently enough to explain the observed amount 
of structure at those scales when the observations are 
interpreted within the framework of a cosmological con- 
stant model (3^ | . 

At present study we are focusing on an alternative 
gravity model featuring deviations from general relativity 
at the scales where observations require to invoke dark 
matter. The question we then ask is not "what grav- 
ity models are consistent with the CDM cosmology" , but 
rather "what properties of GDM are required for cosmo- 
logical viability of a given modified gravity model ?". 

To begin, we will derive an evolution equation for the 
inhomogeneities in a general fluid. To characterize small 
perturbations about the background we can without loss 
of generality adopt the longitudinal Newtonian gauge (36l . 

which is defined by including the two gravitational 
potentials \& and <& to the Robertson- Walker line-element 
as 

ds 2 = a 2 (r)[-(l + 2^)dr 2 + (1 - 2<%g } <fa:W], (2) 
(3) 

where is the spatial 3-metric which in a flat universe 

reduces to the Dirac delta function g9j' — The com- 
ponents of the energy-momentum tensor for a general 
fluid including scalar perturbations can be written as 

n = -p(i+s»), ( 3 ) 

T°i = -(p + pX, (4) 
T) = p(w + c 2 5 N )5) + (Vv, + ^V fc V fe ) IT. (5) 

where an overbar means the background value, 5 N is the 
overdensity and v N the velocity potential evaluated in the 
Newtonian gauge. A nonzero equation of state w = p/p 
would imply that the properties of dark matter differ 
from CDM already in their effects to the overall expan- 
sion of the universe. Then the pressure perturbation 5p 
could differ from zero as well. Our description is not com- 
pletely general as we do not now allow entropic pressure, 



but assume 1 Sp = c 2 Sp, where c 2 = p/p is the sound 
speed (an overdot denoting a derivative with respect to 
the conformal time r). The anisotropic stress is con- 
structed from the scalar potential II with the aid of the 

(3) 

covariant derivatives of the metric g\j ■ The last term 
in Eq.([3]) can appear only at the perturbative level in a 
Fricdmann-Lemaitre universe. In the following, it will be 
useful to employ the gauge-invariant variable 



5 = 5 N + 3H(l + w)^—, (6) 
k 

This 5 is equal to the fractional overdensity in the comov- 
ing gauge, where the velocity potential vanishes. On the 
other hand, it is proportional to the velocity potential in 
the uniform-density gauge [37l ]. 

Following the method of [3l|, one may then derive a 
generalized evolution equation for 5 in the case where 
that energy- momentum tensor is given by Eq.([3]). Here 
we report only the result, which is 

5 = D X H5 + (D 2 H 2 + D k k 2 ) S + P^U + P 2 H 2 U, (7) 

where the dimensionless coefficients are given explicitly 
in the Appendix [AJ Eq. ([7]) determines completely the 
behaviour of cosmological fluctuations with freedom to 
choose any /, u>, c 2 or II. Once the comoving matter 
fluctuation 5 is solved from this second order differential 
equation, the corresponding matter variables in any other 
gauge, as well any metric perturbation, in particular the 
potentials and $ in Eq.0, are uniquely fixed by the 
solution for <5, and can be found as usually by plugging 
this solution into the gauge transformation or constraint 
equations determining the relationships between the dif- 
ferent variables. When w = c 2 = II = 0, Eq. (|7j) reduces 
to the cases considered in [3l| . The term of main inter- 
est to us here is the gradient term Eq. (|A4[1 , whose effects 
easily spoil the agreement of these models with obser- 
vations. This has been extensively discussed elsewhere 
[3ll l32l . I331 . [H, HI] and here we will just discuss the pos- 
sibilities of alleviating the effects of this scale-dependent 
term. Let us though note that its presence stems from 
the new effective matter sources peculiar to the particu- 
lar class of modified gravity models. It has been recently 
found that such a scale-dependent term is absent in vari- 
ous modified gravity theories within the metric formalism 

eh mm m. 



We assume this for simplicity since if to = (which is satisfied 
by baryons and CDM), the entropic pressure should identically 
vanish. On the other hand, a small sound speed described by c 2 
and negligible for the background, could be conceivable for dark 
matter even when w = (this is indeed the standard description 
of baryons). 
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III. A SPECIFIC EXAMPLE ESCAPING THE 
CONSTRAINTS 

One can notice that the dynamics of perturbations in 
a perfect fluid (with II = 0) can be scale-independent in 
the case that the sound speed of the fluid would obey 



ldlogF 
6 dloga ' 



(8) 



where we have defined F = df /dR. In the following we 
will however study an imperfect case with II 0. 

We then use a specific expression for the shear stress 
potential of the cosmic fluid: 



n 



F 



1 - 3K/k 2 \ AHF 



—c 

2 



(9) 



We have kept the prefactor (1 — 3-fT/fc 2 ) -1 here for the 
sake of generality, though it reduces to unity for scales 
smaller than the curvature radius of the universe, and in 
a flat universe it equals identically 1. The first contribu- 
tion in the parenthesis to II is there in order to eliminate 
the effective matter sources due to modified gravity, and 
the second in order to cancel the sources due to possi- 
ble isotropic pressure in any kind of warm dark matter. 
What one then finds with this input is that Eq.© can 
be written in a simpler form 



6 + D 1 H5 + D 2 H 2 5 = 0. 



(10) 



This verifies our claim that the effect of gradient can be 
cancelled by inherent properties of dark matter. 

We can even consider the case when this happens while 
our dark matter fluid at the background level is com- 
pletely pressureless as usual 2 . Setting w = our result 
will further simplify to Eq. (|A7[) in Appendix [XJ Thus 
we have shown that it is possible to avoid the tight 
constraints from structure formation, while keeping the 
background expansion exactly the same as in the CDM 
scenario. In fact, as the evolution in Eqs. (|10IA7p is scale- 
invariant, the shape of the matter power spectrum is ex- 
actly the same for any choice of w, c 2 , and in particular 
for any function f(R). Therefore, if the normalization is 
left arbitrary, no constraints at all arise from comparison 
with the shape of the observed matter power spectrum. 
However, the linear growth rate would be a potentially 
useful test for these models. 

To illustrate this, we will study the toy model / ~ 
R n . This model allows an analytical treatment and its 



2 Though our present approach is phenomenological and not aimed 
at explaining the origin of dark matter stress n, we might note 
that cosmic media satisfying effectively w = but II ^ has 
indeed been considered |43f . (See also |44| for a recent discussion 
anisotropically stressed cosmological fluids and 36] for an exten- 
sive review of structure formation with generalized dark matter). 



predictions have been compared with the data on the late 
time cosmological expansion [2(1 [28l [35l [45l| . In fact, the 
background is simply described by an effective equation 
of state w e ff = — l + (l + iu)/n when K = 0, so that then 



H 2 = 



H 
H 



4n 2 



1 



(3(l + w)-2n) 2 t 2 ' 1 
3(l + w)\ u2 



2n 



1 ( n 3(1 +w) 

2 I 



H 6 



It is also easy to see also that 
F = 3(l + «>)(l-n) 



F = 3 



w + c 2 + (1 + w) 



2n 



- 1 



(11) 
(12) 

(13) 

(14) 
HF. (15) 



Plugging then the expressions (flTj) and (fT4|) in Eq. (fl0|) 
one finds that, keeping w constant for simplicity, 



6 = 



3 (1 + w) - n (1 + 9w) 6 6 (n - 2) (1 + w) 8 



2n-3(l+iu) 



2n-3(l 



(16) 

This admits two power-law solutions, the other one corre- 
sponding to a decaying and the other to a growing mode. 
Setting further w = we find that the latter is charac- 
terized by the growth rate 



_ dlogS \n- |12-7n|| 



dloga 



in 



(17) 



There is an observational estimate [46[ for this r at the 
redhift z = 1/a - 1 = 0.15 from the 2dFGRS data 0, 
implying that r = 0.51 ± 0.11. For the f(R) - R n 
model considered here, an analysis combining the su- 
pernova and baryon oscillation scale constraints shows 
that the background data prefers the values of the ex- 
ponent n = 2.6 ± 0.3 [28[ . This corresponds to about 
0.20 < r < 0.47, thus agreeing with the 2dFGRS con- 
straints. Note that the linear evolution for this same 
model, but without making the assumption Eq.([5]), is 
in gross disagreement with the data[3l|. One may also 
quantify the linear behaviour utilizing the growth index 
7 defined via [48] 



g(a) =6 1 a 



J "-di og a(n m (ay 



(18) 



This single-parameter minimalist characterization of 
modified gravity effects turns out to depend both on the 
matter density and about the exponent n, whereas the re- 
sult for the growth rate r was independent of the amount 
of matter. We find a simple relation, 7 = log r/ log f2 m . 
As the background data constrains quite tightly O m ss 
0.3 [11], we have 0.62 < 7 < 1.34. There are other useful 
probes of modified gravity effects IJj| , but their detailed 
analysis is left for further studies. 
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IV. CONCLUSIONS 

Inspired by the the theoretical developments [25j| . the 
success of the Palatini-/(i?) gravities within the Solar 
system [l9[ and the recent resolutions of their physical 
implications [ll|, we reconsidered the cosmological via- 
bility of these models in view of their predictions for the 
large-scale structure. 

Previous investigations on the subject had established 
that the constraints from matter power spectrum allow 
only models which are practically indistinguishable (by 
e.g. their b ackg round expansion) from the ACDM model 
[3II . I32I l33l l34j In the present paper we examined the 
robustness of those conclusions to the variations of the 
dark matter scenario. We pointed out that the extremely 
tight constraints from structure formation are valid only 
for the /(i?)CDM models, i.e. when the universe energy 



density is assumed to be dominated by pressureless and 
perfect matter. The instabilities occurring then due to 
modified gravity effects, may be less severe or even absent 
in /(i?)GDM scenarios, i.e. when dark matter is allowed 
to have inherent stresses. The cosmological bounds on 
these f(R) models could thereby be drastically loosened. 

This demonstrates that a revision of the nature dark 
energy could also change our view of dark matter. 
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APPENDIX A: EVOLUTION EQUATION 



The equation ((7j) was written as 

5 = DiH5 + (D 2 H 2 + D k k 2 ) 5 + PiHIi + P 2 H 2 Tl. 
Defining F = df/dR, we can then write the dimensionless coefficients as the following: 

-2FH ^FH 2 (1 + 3c 2 - 6w) + F) + 2F 2 H + FF (2H - H 2 (l + 3c 2 - 6w) 



FH 2 [F + 2FH) D 2 = 2FH 



FH [12H 3 (w - c 2 ) + HH(3c 2 + 3w - 2) + Hj + F ( H - H 2 (l - 3w) 
FF ( -2H 2 + HH + 3H 2 (H - 4H 2 )(c 2 - w)^j - 2F 2 H (ff - H 2 (l - 3w) 

(P + 6c 2 FH 



D k 



3 [F + 2FH 



ft-»f£-n. 



and 



2 (fc 2 - 3K) 6F 2 + 9FFH (c 2 - w) - 2F (3F + F{6H + 9(c 2 - w)H 2 - fc 2 )) 



3FH (P + 2F#) k 2 
In the case that Eq.© holds and w = 0, the evolution equation reduces to 



AFH [ FH 2 - + F ) + 3F 2 H + AFF [H - H 



r \\ FH 

2F 2 H 2 IF + 2FH\ { 
2F 2 H (2FH(H - 2HH) + F(2H - 3H 2 )j + 3F 3 H 2 - 3F 2 FH(H - 2H 2 ) 

FF (-4H 2 F + 2HFH - (3F + 2HF)H 2 ) ]s\. 
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